The quantum Prenkel-Kontorova model: a squeezed state approach 
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The squeezed state is used to study the one-dimensional quantum mechanical Frenkel Kontorova 
model. A set of coupled equations for the particle's expectation value and the fluctuations for the 
ground state are derived. It is shown that quantum fluctuations renormalize the standard map to an 
effective sawtooth map. The mechanism underling provides an alternative and simple explanation 
of dynamical localization in quantum chaos. 
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The Frenkel-Kontorova (FK) model describes an 
atomic chain connected by harmonic springs subjected 
to an external sinusoidal potential. This model has 
been widely used to model the crystal dislocations [Q, 
adsorbed epitaxial monolyaers [Q, and incommensurate 
structures The existence of two competing period- 
icities may lead to a rich behavior of the state configu- 
rational properties of the particles Q . Recent years has 
witnessed the application of the FK model to the study of 
transmission in Josphson junction and atomic-scale fric- 
tion -nanoscale tribology, in which the quantum effects 
is very essential @] . 

To get a deep understanding of the nanoscale tribol- 
ogy, it is very necessary to study the quantum FK model. 
However, in contrast to the classical FK model, up to now 
only a few works have been devoted to the effects of the 
quantum fluctuations in the FK model H(j]]. 

Like thermal fluctuations in classical systems at finite 
temperatures, quantum fluctuations play a very impor- 
tant role in quantum systems with finite h. In particu- 
lar, they become crucial and very important at zero tem- 
perature, when thermal fluctuations vanishe. The study 
of quantum fluctuations becomes an important topic in 
quantum phase transitions H and quantum chaos. 

Among many useful tools in study of quantum fluctu- 
ations, the squeezed state, which is a generalization of 
the coherent state, has been proved to be very useful in 
dealing with many-body problems ]9|JTc|| . In this Letter, 
we shall study the effect of quantum fluctuations in the 
one-dimensional FK model by using the squeezed state 
approach. As we shall see later, a set of coupled equa- 
tions for the expectation value and the flucatuation of the 
particle will be derived for the ground state at zero tem- 
perature. We discover analytically how quantum fluctu- 
ations renormalize the external potential, which leads to 
the transition of the standard map in classical FK model 
to the sawtooth map in the quantum FK model. The 
results are found to be in a good agreement with that of 
quantum Monte Carlo (QMC) method. 

The Hamiltonian operator of the one-dimensional stan- 
dard FK model is, 



Pi ,7/- 



Xi - a) 2 - V cos(go^i) 



(1) 



Here m is the mass of particle, 7 the elastic constant 
of the spring, 27r/go the period of external potential. V 
is the strength of the external potential, a the equilib- 
rium distance between two nearst neighbor particles as 
the external potential vanishes. For convenience, we can 
rescale the variables into dimensionless one and obtain a 
new Hamitonian 
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Xi-fif-KcosiXi) 



(2) 



where K = Vq^/j is the rescaled strength of external 
potential. The effective Planck constant h = h-^=, is 
the ratio of the natural quantum energy scale {Tiujq) to 
the natural classical energy scale {^/q 2 )- u 2 =7/771. 

The position and momentum operators for the ith par- 
ticle are written as 



X = vt(<^ + a i)> 

Pi = iJ\{af - at) 



(3) 



Here, af and hi are boson creation and annihilation oper- 
ators which satisfy the canonical commutation relations: 
[&i,a+] = 5ij, [oi,aj] = and [af,af] = 0. 

The squeezed state |$) is defined by the ordinary har- 
monic oscillator displacement operator acting on a 
squeezed vacuum state, 



!<£>(«, /?)> 



10), 



where 



S(a) 
f(f3) 



aipt&j 



(4) 



(5) 



|0) is the vacuum state and a,|0) = 0. S + (a) — 
-S(a),f+(/3) = - f(J3). For simplicity, in what follows 
we will use the abbreviation |<f>) = |$(o!, /?)). 
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It must be noted that if we set (3 = 0, the squeezed 
state is reduced to the coherent state. As we shall see 
later, the coherent state is not able to allow us to study 
the fluctuations. 

Using |$) as a trial wave function for Hamiltonian (|^), 
we can easily find the expectation values of the coordi- 
nate and the momentum operators of the ith particle H , 



(6) 



It is worth noting that the variables Xi and Pj, and 
Gij and Uy form explicitly canonical conjugates JHJ]. To 
find the ground state of the quantum FK model, we shall 
take a variational approach, and these four variables are 
regarded as variational variables. Variation with respect 
to (w.r.t.) Pi immediately yields Pi — 0, and w.r.t. Xi 
yields, 



Xi = ($\Xi\$) = V '!« + ai), 



Xi+i — 2Xi + Xi-i = Ki s'mXi. 



(13) 



P i = ($\P i \*)=-iJ*{c4-a i 



Fluctuations in the coordinate and the momentum are 
given by 



AX? = (<S>\(X l -X l ) 2 \<S>)=hG l 
AP?= ($|(A-A) 2 |$), 

-f— 
T i/i 
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= h(-*- + 4^n 4i G ;fe n fc4 
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The fluctuation covariance between the ith particle and 
the jth particle is 

AXiAXj = ($\(Xi - XiXXj - = HGij, (8) 

where, Gij and Hij are 



where Ki = Kexp^—^Giij, which determines the ex- 
pectation value of the particle's coordinate. Unlike its 
classical counterpart {% — 0, Ki — K), this equation is 
coupled with the quantum fluctuation by hGa. Varia- 
tion w.r.t. liij leads to AhGijUji — 0. Since the fluc- 
tuation Gij cannot be zero, they are always positive, 
we have ILji — 0. To obtain equation for Gy, we first 
(7) take variation w.r.t. Gik and note the following rela- 
tion: = SikSji, where, 5ik and Sji are Dirac delta 
functions. We then multiply both sides of the equation 
by Gkj and take summation over k. Finally we get the 
closed equations for the covariance Gij , 



Gi 



G 



i+ij, 



where 



G^ = i(cosh 2 ^Jp+ + sinh 2 ^Wh + \{M(3 + /3+M) tf , 
Il ij = lGr i 1 (Ml3-f3+M) 
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where, 



M 



sinh VP/3+ cosh y/f3fF 
7W 



(9) 



(10) 



cos X, 



(14) 



(15) 



Since (3 is a symmetric matrix, Gij = Gji and Hij = Hji. 
Furthermore, using the following very important relation, 



($|cosXj|$) = exp --Gu cosX u (11) 



we can finally obtain the expectation value of the Hamil- 
tonian H . 



H= ($|#|$) 



This is a set of equations determining the quantum fluc- 
tuations of the particles, G — {Gij}. G is a N x N 
symmetric matrix which provides all the fluctuation in- 
formation. Its diagonal elements give the variance of each 
particle, while its off-diagonal elements give the covari- 
ance between particles, from which we can calculate the 
correlation function of the quantum fluctuation. These 
equations are coupled with the expectation value Xi. 

Up to this point, we have obtained N x (JV + l)/2 + N 
equations for all variables. These equations provide a 
qualitative picture about the system before we proceede 
to do any detailed numerical analysis. In fact, if we in- 
troduce a new variable, Jj+x = -Xi+i — Xi, Eq.(13) can 
be cast into the map, 



i V ' Lk J 

^(Xi+i-Xi-tf 

i 

2 (^^" Gi+ii+i) — ZfiGi+ii^ 



h+i = h + Ki s'mXi 
Xi+i = ij+i + Xi. 



(16) 



In the same manner, by denoting Qj+xj = Gj+xj — Gij, 
we can also write Eq.(|l4|) into the form of a map, 



Qt+ij — Qij + (G(F — 2))ij, 
Gi+ij = G^ + Qi+ij. 



(17) 



K exp 



(12) The difference between the classical (ft = 0) and the 
quantum FK model from Eq.(Uq) is readily seen. In the 
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classical case the control parameter, namely the ampli- 
tude of the external potential, does not change with the 
position index i. However, in the quantum case, due to 
the quantum fluctuation, the amplitude of the effective 
external potential which acts on the particle changes from 
particle to particle. Because Gu > for any nonzero K, 
Ki < K, which means that the quantum fluctuation re- 
duces the external potential strength acting on the par- 
ticle. Another important difference is that, in the clas- 
sical case, the coordinates of the atoms in the ground 
state are determined by the standard map, whereas in 
the quantum case they are determined by (N + 1) cou- 
pled two-dimension maps. This makes the quantum FK 
model extremly difficult to deal with analytically. 

Before we turn to the numerical calculation, it is worth 
pointing out that in the case of j3 — in Eq.(Q), GV,- = 1/2 
(for all i, j — 1,2, • • • , N), which is the result of the co- 
herent state theory. It is obvious that this cannot be the 
case for a real quantum FK model. So, coherent state is 
not suitable for the study of quantum FK model. 

We now make some comparisons with the quantum 
Monte Carlo (QMC) method. As mentioned before, to 
find the solution from two sets of equations Eqs.(|l^) and 
( |l4| ) is equivalent to find the periodic orbit in a 2(7V + 1)- 
dimensional map Eq.@ and Eq.@. This is still a big 
problem in the nonlinear dynamics to be solved. Never- 
theless, we can make a numerical test of the Eq.([l3]) to 
see whether this equation can give rise to the "sawtooth 
map" Q. 

In Fig. 1 we show the quantum Monte Carlo results 
(left column) and the results calculated from Eq.jl^) 
(right column) by using the QMC's fluctuation data 
Gu in the supercritical regime (K = 5) with % = 0.2 
for an incommensurate ground state. In our quantum 
Monte Carlo computation, as usual, we use the continued 
fraction expansion for the golden mean winding number 
(y/E - l)/2. Thus, we use Q particles which substrated 
into P external potentials with period of 2tt. Periodic 
boundary condition is used: JfC]. Xq + i — Xi + 2nP. 
The winding number is P/Q. The results shown in the 
figure is for P/Q = 34/55. 

By using QMC calculation, we obtained the expecta- 
tion value of the atom's coordinate, from which we can 
construct the so called quantum Hull function (QHF), 
namely Xi (mode 27r) versus the unperturbed ones 
2iriP/Q (mode 2ir), which is shown in top-left in Fig. 
1. The (/-function which defined by |g] 



g i = K- 1 {X i+x -2X i +X i _ 1 ), 



(18) 



is shown in the middle-left of Fig. 1 from the QMC data. 
The quantum fluctuation Gu calculated from QMC is 
shown also in the bottom- left of Fig. 1. 

To compare the squeezed state results with those 
from QMC, we substitute Gu calculated from QMC into 
Eq.(|l3|) and then compute the expectation value of the 



particles' coordinates by using Aubry's gradient method. 
We then construct the quantum Hull function and quan- 
tum ^-function, which is shown in the right column of 
Fig. 1. 
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FIG. 1. Comparision between the quantum Monte Carlo 
(QMC) results (left column) and those from the squeezed 
state theory (SST) Eq.(|l3|) (right column) in the supercrit- 
ical regime where K — 5 with h = 0.2. 



The results from Eq.(13) (right column) agree sur- 
prisingly well with those from QMC for the quantum 
Hull function as well as the ^-function. The most strik- 
ing feature to be noted is the sawtooth shape of the g- 
function in the supercritical regime. This phenomenon 
was first observed by Borgonovi et al || in their QMC 
computation and has been explained as a tunneling effect. 
Later on Berman et al [Jj recovered this phenomenon 
by using a mean field theory including the contribution 
from quasidegenerate states. In the framwork of the 
squeezed state theory, this quantum sawtooth map is just 
a straightforward result of Eq. (|l3|) , which results from 
the quantum fluctuations. 

Our result demonstrates that the squeezed state ap- 
proach indeed captures correctly and nicely the effects of 
quantum fluctuations. 

Finally, we would like to point out that the mecha- 
nism of the reduction of the effective potential due to the 
quantum fluctuations demonstrated above can be applied 
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to explain the quantum suppresion of chaos and relevant 
phenomena such as dynamical localization in quantum 
chaos. The dynamical localization is a well-established 
fact, it was observed by Casati et al |l2| numerically al- 
most 20 years ago and confirmed recently in several dif- 
ferent experiments such as hydrogen atoms in microwave 
fields and so on |b| . Its underlying mechanism is still not 
completely understood. Here, we shall demonstrate that 
by applying the squeezed state approach to the kicked 
rotator, we could obtain a simple and clear picture of the 
dynamical localization. 

Using the squeezed state, we obtain a map like Eq. |L6| ) 
for the expectation value of the angular variable and an- 
gular momentum. But the equation determining Gu is 
different from Eq. (jl4|) , in this case it can be numeri- 
cally calculated. We found that the fluctuation Gu grows 
quadratically with time (kicks), eventually the strength 
of external control parameter Ki becomes very small, 
thus the classical chaos is completely suppressed and 
leads to the dynamical localization. This gives us an al- 
ternative explanation and a very simple picture of the dy- 
namical localization. In turn, it shows that the squeezed 
state is a very useful tool in study the phenomena related 
with the quantum fluctuations. 

In conclusion, we have derived a set of coupled equa- 
tions determining the expectation values of the coordi- 
nate and the quantum fluctuations by using the squeezed 
state as a trial wave function. The results from the 
squeezed-state theory agree with those from the quan- 
tum Monte Carlo method quite well. Furthermore, the 
squeezed-state results give us a very clear understanding 
of the renormalization of the standard map in the clas- 
sical case to the effective sawtooth map in the quantum 
case. Moreover, the squeezed state approach provides an 
alternative and a simple picture of the dynamical local- 
ization observed in many quantum systems. 
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